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Abstract 

We examine the possibility that energetic Standard Model particles contained in the jets produced by 
active galactic nuclei (AGN) may scatter off of the dark matter halo which is expected to surround the AGN. 
In particular, if there are nearby states in the dark sector which can appear resonantly in the scattering, 
the cross section can be enhanced and a distinctive edge feature in the energy spectrum may appear. We 
examine bounds on supersymmetric models which may be obtained from the Fermi Gamma-ray Space 
Telescope observation of the nearby AGN Centaurus A. 



I. INTRODUCTION 



Active galactic nuclei (AGN) are among the most powerful natural accelerators in the Universe, 
producing jets of charged particles which travel macroscopic distances and can even be imaged 
from the Earth through the subsequent production of high energy photons. While the nature of 
the jet constituents (proton or electron) remains a subject of some controversy, there is no doubt 
that they are among the highest energy and most interesting objects in the cosmos, with energies 
reaching into tens of GeV and with huge luminosities. 

The supermassive black holes which power these jets are typically hosted in very large galaxies, 
where one would expect the largest concentrations of dark matter (DM) also occur. Taken together, 
the combination of dense dark matter with high energy particles represent a unique opportunity to 
study the interactions of dark matter and Standard Model (SM) particles at high energies. Indeed, 
while particle accelerators on the Earth reach higher energies and with (much) better understood 
beams, the huge density of dark matter expected close to an AGN nevertheless renders it a unique 
opportunity to learn about dark matter interactions with high energy SM particles. 

If there is some other particle in the "dark sector" which can be produced as a resonance in the 
scattering of dark matter with a jet particle, the rate for scattering in the high energy tail of the 
jet may be enhanced. Beam particles with energies E greater than, 

9 o 

1 m| — m~ 

E>-^ * (1) 

~ 2 m x 

(where m^- is the mass of the dark matter and mg the mass of the resonance) have sufficient energy 
to produce the resonance on-shell, leading to an enhancement of the cross section. Since the beam 
particles are presumably charged (so that they can be efficiently accelerated by the AGN), the 
resonant state should also be charged. As a result, its decays can produce gamma rays which 
can stick out from the background produced by the AGN and jet themselves. Such a situation is 
common in the most popular theories of dark matter, such as in supersymmetric theories in which 
a neutralino plays the role of dark matter pQ , and which can be up-scattered by an energetic beam 
electron into a selectrorj^] Such a degeneracy may even be central to realizing a neutralino 
WIMP as dark matter, to allow for efficient enough coannihilation to produce the observed density 
from thermal freeze-out. In fact, a very degenerate selectron is challenging to detect at the LHC, 
because its visible decay products have energies controlled by nig — rriy, and are typically very 



1 Similarly, in theories with Universal Extra Dimensions [3], the dark matter particle may be a Kaluza-Klein 
excitation of the hypercharge vector boson [HE], which may up-scatter into a Kaluza-Klein mode of the electron 
0. 
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soft. Thus, the signal from AGNs is enhanced precisely where the LHC has difficulty, and can help 
study motivated but otherwise experimentally difficult regions of parameter space. 

Since the resonant state must be charged, its subsequent de-excitation includes higher order 
processes in which a photon is produced in addition to the primary SM and dark matter decay 
products. While suppressed by a, such processes may be enhanced by collinear logarithms that 
rise when the photon is emitted along the same direction as the final state charged SM particle. 
Unlike charged particles, these radiated photons have very long mean free paths, and could be 
detected by detectors on the Earth as an indicator that the jet-DM scatterings are taking place. 
For dark matter masses on the order of a few hundred GeV, these photons fall precisely in the 
range for which the Fermi Gamma-ray Space Telescope is sensitive, and thus it is precisely the 
right time to explore such signals. In fact, as pointed out in Ref. jB], the rise in cross section as 
beam particles cross the e threshold leads to a very interesting feature in the spectrum of gamma 
rays, a rise in the number of high energy photons followed by a sharp edge. This feature provides 
an interesting discriminant from the background of gamma rays emitted from the jet itself, one 
which conventional astrophysical sources are expected to have difficulty faking. 

This article is organized as follows. In Section [XTJ, we discuss the astrophysical inputs to the 
computation, namely the density profile of dark matter around the AGN and the intensity and 



energy spectrum of the beam particles. In Section III, we revisit the computation of the cross 



sections for scattering in supersymmetric theories. In Section IV we examine the bounds one can 



place for different astrophysical assumptions. We conclude in Section |V} 



II. ASTROPHYSICAL INPUTS 



In this section, we discuss the astrophysical inputs to our estimates. These boil down into the 
properties of particles inside the AGN jet, and the density of dark matter surrounding it. In both 
cases, we follow closely the discussion of Ref. [6], which contains a more detailed discussion, and 
study of the impact varying these parameters has on the signal. We restrict ourselves to looking 
at Centaurus A, which offers a good balance of being relatively close by (D ~ 3.7 Mpc |7J) as well 
as a massive object surrounded by a dense halo of dark matter. 
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A. Dark Matter Distribution 



The impact of the dark matter distribution is through its integral along the "line of sight" of 
the jet particles, 

r''max 

5dm = / dr p(r) (2) 



where p(r) is the density of dark matter at a radial distance r from the center of the AGN, r m ; n 
characterizes the minimum distance of interest (the base of the jet), and r max the distance at which 
the jet loses cohesion and becomes irrelevant. In practice, results depend sensitively on r m ; n and 
are rather insensitive to r max [6]. 

The dark matter density distribution surrounding a compact object was modeled in the colli- 
sionless limit and under the assumption that the central black hole grew adiabatically by accretion 
in Ref. [8]. The resulting steady-state distribution is characterized by a central spike of dark matter, 
which gets depleted in the core by pair annihilation. The resulting distribution is thus sensitive 
to the underlying particle physics model of dark matter, with the central density enhanced for 
WIMPs with very low (or in the case of a purely asymmetric dark matter halo, zero) annihilation 
cross sections. Ref [8] begins with a cusped density profile pi{r) = po(r /ro)~ 7 , and determines a 
final profile (close to the center) of the form, 

P[r) = P '(r) + Pcoie ' (3) 
where the core density is given in terms of the velocity-averaged annihilation rate (av) and the age 
of the black hole £bh as, 

(4) 
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The dependence on r is contained in the function, 



where 7 sp is slope of the density profile in the dark matter spike, R sv its radius, and Rs is the 
Schwarzschild radius. Below 4i?5, the dark matter is accreted into the black hole, and has p = 0. 
The spike radius and power are given by, 

R sp = « 7 r y-^J and 7sp = , (6) 

where a 7 ~ 0.2937 4 / 9 for 7 <C 1 [8]. This analysis neglects the influence of stellar populations 
close to the the supermassive black hole, which can gravitationally scatter dark matter out of the 
spike, resulting in a somewhat shallower profile [9j. 



Following Ref . [6] , we normalize the distribution such that it falls within the uncertainty on the 
measured mass of the black hole itself. The mass of the black hole is determined by measurements 
of stellar kinematics to be |1U| . 

M BH = (5.5 ± 3.0) x 1O 7 M (7) 
which corresponds to Rs — 5 x 10~ 6 pc. Given a choice of 7, we determine po by requiring, 

rl0 5 R s 

3 x 1O 7 M > / dr 47T7- 2 p{r) , (8) 

where W 5 Rs characterizes the distances at which the stellar populations are measured. 

Ultimately, the factor 8bm is very sensitive to the product (av) x t^H, which determines the 
density within the spike. For black hole ages between 10 8 to 10 10 years and cross sections ranging 
from 10~ 26 cm 3 s _1 (appropriate for a thermal relic) to 10 -30 cm 3 s -1 (typical for the "coannihi- 
lation" region of the MSSM, in which the superpartners of the charged leptons play an active role 
in freeze-out), one finds <5dm ranging from, 

5 DM ~ 10 8 to 10 11 M Q I pc 3 (9) 

with rather mild dependence on the choices of 7, r m ; n , and ro ~ r max ~ 15 kpc [6]. 

B. AGN Jet Particles 

For our purposes, the details of the AGN jet geometry are not very important, while modeling 

the energy distribution of particles in the jet is crucial. We assume most of the particles in the 

Centaurus A jet are electrons. Based on observation of emitted gamma rays by the Fermi LAT 

experiment we assume the distribution in the electron boost is a broken power law, 

-1 

(Tmin ^ 7 ^ Tmax ) • (10) 

Here the primed variables refer to the "blob frame" , in which the electrons move isotropically. The 
parameters s\, s 2 , j' br , and j' max are fixed to [TT] 

Sl = 1.8, S2 = 3.5, 7 ^ = 4xl0 5 , 7min = 8xi0 2 , 7^ = 10 8 . (11) 

The normalization k e can be determined from the jet power in electrons, which is defined in the 
black hole frame as, 

L * ^ L Trt^M) L d7(me7) "ST ( 7(rs(1 " M) ) ' (12) 

(13) 
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where /U = cos 9 (9 is the polar angle with respect to the jet axis) and 7 = e/m e where e is the energy 
of the electron. Similarly, the primed variables in the blob frame are // = cos#' and 7' = e'/m e . 



The quantities in the two frames are related by the blob velocity /3g and boost Fb = — (3% 

( ~ 3 for Centaurus A) as, 



B 



or equivalently, 



I- Pb H 



7' = 7 Tb (1 - /?£}//) , 
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The limits on the 7 integral are given by, 
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(14) 



(15) 
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whereas the limits on /x in the black hole frame correspond to 0.9 — 1, leading to a highly collimated 
jet of electrons. 



The function d&i AGN ^ /d*y can be written explicitly as, 



d$ (AGN) / x 1 

^ ^(r fl (i-/3 flA1 ))J =-fc e 



7r B (i -Pbij) 



1+ 7(r B (l-/3 B M))/7fcr- 



Saturating the upper limit of L e < 10 erg s 1 [BJ, we arrive at 

1 



-k P = 4.2088 x 10 4y s 
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(17) 



(18) 



III. CROSS SECTIONS 



The total cross section for the process e x ~ > e X 7 may be written, 

*r = A- A - r)^^^^ZlM? , (19) 

where the \ Ai\ 2 is the sum of all matrix element squared, including both s- and u-channel 
exchange of e, shown in the Fig. [TJ We ignore i-channel contributions mediated by the SM Z 
boson as it is non-resonant, and thus negligible. The amplitudes corresponding to these diagrams 
are labelled Mi - M 6 . 

The experimental observable is the distribution of photon energies, vS v = E^^jr- at a fixed 
angle relative to the jet axis (~ 68° for the Centaurus A). For fixed incoming electron energy, we 
integrate over the phase space of the final state neutralino and electron, resulting in d 2 a / dE^dQ^. 
We work in the rest frame of the initial neutralino dark matter particle, which is expected to be 
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FIG. 1: Feynman diagrams for the process e x — s- e x 7 through e exchange. The amplitudes 
corresponding to these diagrams are, sequentially, Mi — Mq. 

moving non-relativistically (~ 300 km/s [12J) in the rest frame of the Earth, d 2 o / dE^dQ,^ can be 
expressed, 



d 2 a 1 1 

dE^dVL^ ~ (2vr) 2 4m^E 2 8E' 



(20) 



where J is a Jacobian factor from the integration of the function <5 4 (pi +P2 —p'± —P2~ r ) over dE' 2 
and can be explicitly written as 

E' 2 (pi + p 2 - r) ■ 2 



J 



E i E[E' 2 
E' 2 (P2x ~ r x )p' 2x + (P2y - r y )p' 2y + (p 2z - r z )p' 2z 



E[E' 2 



(21) 
(22) 



The full expression for the amplitude squared may be found in Appendix [A] The cross section 
is enhanced in three kinematic configurations: 

• Mi_3 are resonantly enhanced when the incoming electron energy is such that the interme- 
diate selectron is approximately on-shell [2|. 

• M 2 and M5 have a collinear enhancement when the emitted photon lines up with the final 
state electron [6]. 



M 2 and M5 have a soft enhancement for low energy photon emission. 
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Note that because the jet axis makes a fixed angle with respect to the line of sight from the Earth, 
there is no possibility of a collinear enhancement from the initial electron in M\ and M4. In practice, 
the soft enhancement (the infinity of which is formally matched by the one-loop QED correction 
to e~ x ~ > e ~ X scattering) is not observationally interesting. The dominant contribution comes 
from the on-shell selectron, collinear photon region of M 2 : 

2 _ a cS (pi • r)(pi -p 2 ) 



\M 2 \ 2 = ^ ctt V 1 — \ , (23) 

2 



|S S | 2 l(p' 2 -r] 



where a cB = 4ira EM (af j + a^) is the effective neutralino-selectron-electron coupling. We assume the 
left- and right-handed selectrons are degenerate in mass, and a mostly bino neutralino, for which, 

4 = 2 5 2 tan0 2 , (24) 
a 2 L = a 2 J A . (25) 

The factor 1/|X S | 2 is the selectron propagator, which has been modified to maintain gauge invari- 
anc^J using the prescription of Ref. [14] . To avoid the divergence from soft radiation, (p' 2 • r) in 
the propagator has been shifted by m 2 e /2. 

The differential cross sections d 2 a j dE e d£l^ are shown as a function of the incoming electron 
energy in Fig.[2j for fixed photon energies of 10 GeV, 30 GeV and 50 GeV. As a benchmark point, 
we have chosen the mass of the neutralino to be = 60 GeV and the mass of the selectron to be 
rrig = 100 GeV. Notice that for a photon energy of 10 GeV or 30 GeV, there are two peaks which 
arise from the (shifted) resonant behavior of M 2 and Mi separately. However, as the photon energy 
increases, the incoming electron energy must be large enough to produce such a hard photon in the 
scattering process, closing off the possibility for scattering (for fixed E^) below a certain threshold. 
Once the threshold rises above ~ Mg/2, it becomes impossible to produce a final state photon from 
an on-shell e decay, and the M 2 resonance effectively disappears. 

IV. BOUNDS ON SUPERSYMMETRIC PARAMETER SPACE 

Putting the astrophysical inputs together with the cross section, we have: 

— = [5 DM ) x J dE e [ W2 ^) x {-J^y (26) 



II A 



II B 



where the three factors 5dm, d^ GN / dE e and d a/dE^dQ,^ have been discussed in Sections 
and III, respectively. We adopt two choices of 5dm = 10 10 M pc~ 2 and 5dm = 10 11 M pc -2 and 



Proof of gauge invariance can be found in Appendix B 
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FIG. 2: Differential cross sections, d 2 a/dE~ / dQ~ / for the process e x — >• e X7, for a neutralino 
mass of 60 GeV and selectron mass of 100 GeV, and fixed photon energies (top to bottom) of 

10 GeV, 30 GeV, and 50 GeV. 



to begin with, the benchmark point with a bino neutralino of mass rriy = 60 GeV and degenerate 
left- and right-handed sleptons of masses mg = 100 GeV. In Figure [3j we plot the distribution of 
photon energies for the benchmark point, for both values of <5dm discussed above. We note the 
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vS v 




E y (GeV) 

FIG. 3: The photon energy distributions for (5dm = 10 10 Mq pc -2 (lower curve) and 
<5dm = 10 11 Mq pc~ 2 (upper curve) exhibiting a drop-off around E 1 ~ 35 GeV. 

drop at energies around Ely ~ 35 GeV, where, as discussed above, the required photon energy is so 
large that it is no longer produced efficiently from on-shell intermediate selectrons. 

Using the existing energy sprectrum data from Fermi LAT observations of Centaurus A, we can 
put bounds on the coannihilation region of SUSY parameter space. Ideally, one could do a shape- 
based analysis, as advocated in Ref. |6j, looking for a significant drop in the energy spectrum, which 
occurs as one demands the photon energy be high enough that kinematics forbid it from arising 
from an on-shell selectron decay. However, with the limited statistics of the current data, it is more 
feasible to consider a counting experiment analysis in each energy bin. We place 95% confidence 
limits on regions of MSSM parameter space where a contribution larger than the expected 95% 
fluctuations (based on the error bars) in the spectrum is predicted. We use the quoted uncertainties 
on the Fermi measurement points, 

VB ~ 2.5 x 1(T 12 erg s" 1 cm" 2 , (27) 

roughly independent of energy [TT]. These bounds are subject to our assumptions concerning 5dm 
and the beam composition and energy spectrum. 
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FIG. 4: Contours of 95% CL constraint on the — Am plane based on one year of observation 
of Centaurus A by the Fermi LAT. The colored regions indicate constraints from each energy bin 
reported by Fermi, and the black dashed line indicates mg = 100 GeV. The two panels assume 

different 5dm, as indicated. 

In Figure [4], we plot the 95% CL bounds in the plane of and Am, where 

Am = mg — m^. (28) 

is the mass splitting between the neutralino and selectron. The upper panel corresponds to <5dm = 
10 10 Mq pc -2 , and the lower to 5dm = 10 11 Mq pc~ 2 . These values correspond to a neutralino 
annihilation cross section of (av) ~ 10 -30 cm 3 s _1 (appropriate for a coannihilation scenario) and a 
black hole lifetime of £bh ~ 10 8 or 10 10 years, respectively. We continue to assume a bino neutralino 
and selectrons which are degenerate in mass. Since the optimal energy bin for the search varies 
with m^ and mg, we derive the constrained region from each energy bin independently, as indicated 
on the figure. We find that the largest resolving power comes from the highest energy bins, but 
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nevertheless the lower energy bins provide interesting constraints, and a combined shape analysis 
would probably do a little better than treating each bin independently. The black dashed line on 
each panel indicates the line of mg > 100 GeV, roughly the bound from LEP-II null searches [T5]. 

V. CONCLUSIONS AND OUTLOOK 

Active galactic nuclei are among the most energetic natural accelerators in the Universe. They 
are also located in regions rich with dark matter, and provide a unique opportunity to study the 
high energy interactions of WIMPs with SM particles. In this paper, we have estimated the regions 
of MSSM parameter space which can be probed from Fermi observations. We find that there is 
greatest sensitivity for models with low annihilation cross sections and nearly degenerate sleptons 
- as occurs in coannihilation regions of the MSSM, where sleptons are active during freeze-out. 
Indeed, sleptons are very challenging to discover at the LHC under ideal circumstances, and almost 
impossible when nearly degenerate with the lightest supersymmetry particle. Observations of the 
gamma rays from AGNs can help cover an interesting and important gap in LHC coverage of 
MSSM parameter space. 

Ultimately, the limiting factors in utilizing these naturally occurring accelerators as probes 
of the property of dark matter are the uncertainties in the properties of the AGNs themselves 
and the dark matter environment around them. These uncertainties can be mitigated by further 
observations, to help pin down the underlying properties of the jet, and observation of the dynamics 
of the surrounding stars and gas to better constrain the gravitational dynamics. The MSSM signal 
considered here also leads to a distinctive feature in the gamma ray spectrum, which could be 
exploited with more statistics to minimize the dependence on the unknown background emission. 
Ultimately, if the particle properties of dark matter can be pinned down more precisely using data 
from the LHC or direct or indirect dark matter detection experiments, one could even turn the 
process around and use collisions between jets and the WIMPs to learn about the astrophysics of 
these fascinating systems. 

VI. ACKNOWLEDGEMENTS 

T. Tait is glad to acknowledge conversations with L. Costamente, S. Profumo, and L. Ubaldi, 
and the hospitality of the SLAC theory group, for their generosity during his many visits. The 
work of AR and JH is supported in part by NSF grants PHY-0653656 and PHY-0709742. The 



12 



work of TMPT is supported in part by NSF grant PHY-0970171. 



Appendix A: Matrix Elements for e \ ~~ ^ e XI 

In this appendix, we summarize the full expression for the amplitude squared for the process 
e~x ~~ ^ e ~X"f- We begin with some notations, 
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Here pi , ,p2 ,p' 2 ,r are the four- vectors for the momenta of the incoming and outgoing dark 
matter particles (neutralino), incoming and outgoing electrons and photon. Ei s , i? 2s , Ei s i, and 
E 2s i are the energies of the initial and final neutralino and electron in the lab frame where the initial 
neutralino is at rest. The s and s' subscripts correspond to the s and s' diagrams, respectively. Tj 



are the momentum dependent decay widths and E s , £ s ' and £ u , S u / are the propagators in the s- 
and u-channels. ^3 is the modified vertex factor in amplitude M3 which insures gauge invariance, 
see appendix [Bj The squared amplitudes, including the interference terms, are: 
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- 2(pi • p 2 ) 2 (pi • p 2 ) - (pi • P 2 )(P2 • P 2 )(p 2 • r) - (pi • p' 2 )(p 2 • r)' 
+ (Pi • P 2 )(p 2 • r)(pi ■ p 2 ) - (pi • p' 2 )(p 2 • r)m| 



-2a cft - 



sp^a^wg72) 1 (pi " p ' 2 ^( p2 " P2) + (pi " p ' 2 ^( p2 " r) 

- (Pi • P 2 ) 2 (Pi • P2) - (Pi • P 2 )(p 2 • P 2 f + (Pi • P 2 )(p 2 • P 2 )(Pi • r) 

- 2(pi • p' 2 )(p 2 • p' 2 )(p 2 • r) + 2(pi • p' 2 )(p 2 ■ p' 2 )(pi • p 2 ) - (pi • p' 2 )(p 2 • p' 2 )m\ 
+ (Pi • P 2 )(Pi • r)(p 2 • r) - (pi ■ p' 2 )(pi • r)(pi • p 2 ) - (pi ■ p' 2 )(p 2 ■ r) 2 

+ 2(pi • p' 2 )(p 2 • r)(pi ■ p 2 ) - (pi ■ p' 2 )(p2 ■ r)m\ - (pi ■ p 2 )(pi ■ P2) 2 
+ (Pi • P 2 )(Pi • P2)m\ - (p 2 ■ p' 2 ) 2 m\ - (p 2 ■ p 2 )(pi • r)m\ 

- 2(p 2 • p' 2 )(p 2 • r)m\ + 2(p 2 • p 2 )(pi ■ p 2 )m| - (pi ■ r)(p 2 ■ r)m\ 
+ (Pi • r)(pi -p 2 )m\ - (p 2 ■ r) 2 m 2 - + 2(p 2 ■ r)(pi ■ p 2 )m\ - (p x ■ p 2 ) 2 m\ 



(A3) 
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2Re(M 1 M{ 



2Re(M 1 Ml) 
2Re(M l M\) 



2Re(M A M\) 
2Re{M 2 M\) 
2Re{M 2 M\) 



2Re(M 5 Ml) 



Re 
Re 
Re 



S s ,S„fe. r ) [(Pi ■ P 2 ) + (P2 • P 2 ) + (Pi • r ) + (P2 • r) ~ (pi ■ p 2 ) 
S s ,S M ,( P2 .r) C ((^.r)-^/2) [fa " I^)(P2 " P 2 ) + " P 2 )(pi " r) ~ ( P 2 ■ P 2 )(p 

2(pi • p 2 )(p 2 • p' 2 ) + (Pi • P 2 )(P2 • r) 



P2) 



S s ,S u ,S u (p 2 -r) 

- (Pi • P 2 )(Pi • P 2 ) + (p 2 • P 2 )(Pi • r) + (P2 • P 2 )(p 2 • r) 

~ (P2 • P 2 )(pi • P2) - (pi • r)(pi • p 2 ) - (p2 • r)(pi • p 2 ) + (pi • p 2 ) 2 



Re 
Re 



— Q; 772 

S B S;( P2 ,)((^.r)- m ^/2) (fa • ?4)(P2 • P 2 ) + (P2 • P 2 )(Pl • r) - (pa ■ P 2 )(pi ■ P 2 ) 

[a cff m|(p2 • r)]/[E a £ u /((p£, • r) - m*/2)] 



<*eff m- 



2i?e(M 4 M 3 t ) = i?e 



^(Pi • P2XP2 • P 2 ) + (Pi • P 2 )(P2 • r) 

" (Pi • P 2 )(Pl • P 2 ) - (P2 • P2) 2 ~ (P2 • P 2 )(P 2 • 

- (pi • r)fa - P2) - (P2 • r)(pi -p 2 ) + (pi ■P2) 2 

"off m| h 



%^r) ((Pi ■ P 2 )(P2 • r ) - (Pi • P 2 )(Pl • P2) 
" (P2 • P 2 )(Pi • r) - (P2 • P 2 )(p 2 • r) + (P2 • P 2 )(Pl • P 2 ) 
- (pi • r)fa - P2) - (P2 • r)(pi -p 2 ) + (pi ■P2) 2 
Re 



-a cB mi f 3 



S„,S;S;,((^-r)-ml/2) ^0»1 • ^)(P2 " ^) + (Pi " P 2 )fa " 
(Pi • P 2 )fa • P 2 ) + (P2 • P2) 2 + (P2 • P 2 )(P 2 • 

(pi • r)fa - P2) - (P2 • r)(pi -p 2 ) + (pi -p 2 ) 2 



2Re(M 3 Ml) = Re 



-2a cB ml f* 



(pi -p 2 )(p 2 -p' 2 ) + (p 2 -p 2 ) 



/ \2 



+ (p 2 • P 2 )(Pi • r) - (p2 • P 2 )(Pi • P 2 ) - 2(p 2 • p 2 W 



(A4) 



Appendix B: Proof of Gauge Invariance 



In this appendix, we verify the invariance of the amplitude under U (1)em gauge transformations. 
We begin with the set of u-channel Feynman graphs. The amplitudes of M^M^ and M§ can be 
written as, 



M 4 



x(Pi)(a L P L + a R P R ) 



7- fa -r)-m, 



-{-ieY)e(P2) 



e(p' 2 )(a L P L + a R P R )x(pi) 
x{p'i){a L P L + a R P R )e{p 2 ) 



(p 2 — Pi — f) 2 — m, 

i 



(p 2 - Pi) 2 - m| 
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q - pi - P2 




e(vo) 



e(v'n) 



FIG. 5: A general "bubble" diagram containing no photon radiation. 



7'(P2 + r ) — m e 



(a L P L + a R P R )x(pi) 



M 6 



x(v'l)(a L P L + a R P R )e(p 2 ) 



e(p' 2 )(a L P L + a R P R )xipi 



;(-ie)(2p 2 -2p' 1 -rY- 



. (P2 - K) 2 ~ m l ' (P2 — Pi - "r) 2 - 

Applying the Ward identity to M^-Mq individually gives rise to 



M 4 = e 
M 5 = - 
M fi = e 



x{p'i)(aLPL + a R P R )e(p 2 

x(Pi)(o>lPl + a R P R )e(p 2 
x(p'i)(cilPl + a R P R )e(p 2 



e(p' 2 ){aLPL + a R P R )x{pi) 
e(p' 2 ){a L P L + a R P R )x(pi 
e(p' 2 )(aLPL + a R P R )x(pi 



(p2 - p[ - r) 2 - m 

i 



(p 2 - Pi) 2 ~ m, 



(Bl) 



(B2) 



. (P2 -p'i) 2 - m?~ (p 2 - p[ — r) 2 - m\ 
and they sum to zero as expected. 

The set of s-channel diagrams is somewhat more complicated. To begin with, we consider the 
process with no final state photon, but a sum of bubble diagrams as shown in Fig. [5| which will be 
useful later. The key is the imaginary part of the loop which produces the width of the selectron. 
Writting the sum of diagrams with all numbers of bubble insertions as, 



p 2 — m~ p 2 — m| 



(-i£(p)) 



p2 _ m | 



+ 



p 2 — m 2 — 



we have 



TV 



{o-lPl + a R P R ) — r {a L P L + a R P R ) 

J-q-rrie 7-{<l ~ P) ~ m x 



(B3) 



(B4) 



where the width is determined by mF = Im(S(p)) at p 2 = m 2 . The diagram with no photon 
radiation is further simplified to, 



x(Pl) (a L P L + a R P R ) e(p 2 



p 2 — m~ — irrieT 



e (P2) (o-lPl + a R P R ) xip'x) 



(B5) 
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with p = pi + P2- 

When the photon line is added, we get new diagrams by attaching the photon in all possible 
ways. The first possibility is 




Mi 



x(pi)(a L P L + a R P R )- 



j.(p 2 - r) - m, 
e(P2)(a L PL + a R P R )x(p[ 



\-ie)^e(p 2 ) 



(Pi + P2) 2 - m i - imeT(pi + p 2 - r) 

(B6) 



We have taken the momentum of the photon (with polarization fj,) to be r going out from the 
diagram. This graph corresponds to M\. We also have the diagram corresponding to M2, 






Pi +P2 



q - pi - p 2 




Mo 



x(pi)(cllPl + a R P R )e(p 2 



(pi + P2) 2 - ml- imeT(pi + p 2 ) 
(a L P L + a R P R )x{pi) 



1-{P2 + r ) 

While for M3, the photon attaches to the selectron; there are two diagrams where one is 



(B7) 




Xip'x) 




7(r) 



Pi + P2 - r 



q-pi-pi 




e(v2) 



e(v'n) 
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M 3 , 



x{p\){a L P L + a R P R )e(p 2 ) 
i 



e(p' 2 )(a L P L + a R P R )x(p'i] 



(-ie)(2p2 + 2pi -r) 



. (pi + P2) 2 - m§ - imeT(pi + p 2 ) (P2 + P2 - r) 2 — m| — imgr(pi + p 2 - r) 
and the other one has the photon attached to an internal loop, 



(B8) 




P1+P2- r 




e(w>) 



e(v'n) 



M: 



3.2 



x{P2){aLPL + a R P R )e(p 1 ] 
i 



e{p[)(a L P L + a R P R )x(P2, 



-&{p x +p 2 ,-r) 



. (Pi + P2) 2 —m~ — inrieT(pi + p 2 ) (pi + P2 ~ r) 2 — m| — rni g r(pi + p 2 - r) 
Now we need to extract the imaginary part of — —r) in this vertex, 



i£^(p,-r) = J 



(27T)* 



TV 



r )\q — r)—m e ^y.q — m. 



(a L P L + a R P R ) -, r (a L P L + a R P R ) 

J.[q-p) - 



TV 



(2vr) 4 
(a L P L + a R P R ) 



-(-ie)-y.(-r)- 



< ) I t^4tTV 



j.(q - r) -m e ' ~, • - in, 

i 

l-{q-p) - 

i 



(a L P L + a R P R ) 



(2tt)4 

(a L P L + a R P R ) — . 

i\q-p) - ™>x 

(— e)(— «'E(p - r) + iS(p)) . 



-f.(q-r)-m e ^.q-m e 



{a L P L + o^-Pr) 



(B9) 



(BIO) 



We can write generally, —iT,^(p, —r) = A {p + p')^ 1 + B . The equation above tells us that 

A(p + p').(-r) = (- e )(-ffi(p - r) + »E(p)U = (~ e )(-^(P ~ r ) + ^(p)) _ 

— 2p • r 

The coefficient B is irrelevant since r • e = 0, and the vertex is effectively — — r) = (— ie)(p + 
p' ) m ( ~ s ^J 2 ^ ^ s (p) ) which can be thought of as a correction to the tree level vertex. Since we are 
only interested in the imaginary piece of the correction, 



-ie)(p + p 



r^ \-T(p-r) + T(p)) 



-2p ■ r 



(B12) 
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we can combine the last two diagrams, arriving at the expression, 



M 3 



xipi){a L PL + a R P R )e(p 2 ) 



{pi + P2) 2 -m\- irrieT(pi + p 2 ) 



e(p' 2 )(a L P L + a R P R )x{p'i) 
-ie)(pi +p 2 -r) 



1 + imz 



+ a R P R )e(p 2 ) 



TQi +P2) - TQi +p 2 - r) ~ 

2(pi +P2) ■ (~r) 
i 

(pi +P2 - r) 2 — m| — irrieT(pi +p 2 — r) 

1 



e{p' 2 ){a L P L + ai?-Pi?)x(p'i)(-ie)(pi +p 2 - r)** 



1 



1 



2(pi+p 2 )-{-r)l 
(B13) 



. (pi + P2) 2 ~ m \- imeT(pi + p 2 ) (pi +P2 - r) 2 —rri~ — zm g r(pi + p 2 - r 
Therefore, M\ + M 2 + M3 = 0, as expected, and we complete the proof of the independent gauge 
invariance for both the s- and the u-channel sets of Feynman graphs. 
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